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ABSTRACT 



Aims. We report on an estimate of the bar pattern speed Clp for the low surface brightness spiral galaxy UGC 628. 

Methods. We applied the Tremaine- Weinberg method to high resolution Ho- velocity and integrated emission maps of this dark matter 

dominated galaxy. Observations were made at the CFHT using the optical Fabry-Perot interferometer, FaNTOmM. 

Results. The Tremaine- Weinberg method estimates a bar pattern speed of (11.3 ± 2.0) km s"' kpc"' for UGC 628, which is among 

the lowest values found for a spiral galaxy. The corotation radius R,.. of the bar and the gaseous disc is Re = 9.8;^2"g kpc, implying a 

ratio "R = Rc/as of 2.0^0 3, where as is the bar radius. The ratio is well beyond the usual range of values, l.O <'R < 1.4, found for fast 

bars of high surface brightness barred galaxies. It implies that the bar in UGC 628 is slow. 

Conclusions. As shown through the use of numerical simulations, fast bars survive when the inner mass distribution of galaxies is 
dominated by the baryons over the dark matter Our result suggests that the presence of slow bars in galaxies is likely related to the 
dominance of dark matter over the mass distribution. 
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1. Introduction 

In the last decade, estimates of the pattern speed of stellar bars 
in spiral galaxies have shown that they have high amplitudes 
(Gerssen, Kuijken & Memfield 1999; Aguerri, Debattista & 

■ Corsini 2003; Corsini et al. 2007). As a consequence, the ra- 

] tio of the corotation radius, R^,, to the semi-major axis of the 
bar, aB, is usually observed in the range [1; 1.4], within the un- 
certainties. A bar whose ratio falls within this range is called a 

I fast bar because its pattern speed cannot be larger, according to 
theoretical models. A bar that would rotate faster would have 
a smaller corotation radius at fixed bar semi-length (fi < 1) 
and should not exist. Indeed, analytic solutions of the equations 
of motion and numerical simulations have shown that periodic 

' orbits of stars and gas are aligned perpendicular to the bar in 
regions beyond corotation and tend to destroy it (Contopoulos 
& Papayannopoulos 1980; Athanassoula 1992). The numerical 
models of Athanassoula (1992) also show simulated bars with 
1 .0 < !R < 1 .4, which result is in excellent agreement with the 
observations. Conversely, stellar bars having H > 1.4 are called 
slow bars and can exist, though they have not been observed yet. 

Classifying a bar as fast or slow thus relies on the measure- 
ments of its pattern speed, its extent and on the details of the 
disc rotation curve. Most of pattern speed determinations were 
made through the application of the Tremaine- Weinberg method 
(TW; Tremaine & Weinberg 1984). Starting from the continuity 
equation, Tremaine & Weinberg showed that a pattern speed can 

^ Svdx/ Sxdx where / the disc 
inclination, S is the surface density of the tracer, v is the line-of- 
sight velocity and x the position along slices parallel to the disc 
major axis. Since the continuity equation underlies this estimate, 
applying the TW method to the old stellar population is most ac- 



curate as it is expected to satisfy the continuity equation better 
than the gas. For this reason, authors have generally used long- 
slit spectroscopy of stellar absorption lines of early type barred 
discs for these estimates (Merrifield & Kuijken 1995; Aguerri 
et al. 2003; Gerssen, Kuijken & Memfield 1999, 2003; Corsini, 
Debattista & Aguerri 2003; Corsini et al. 2007). Nonetheless, the 
TW method has been applied with some success to the kinemat- 
ics of the neutral, molecular, and warm ionized gas (Bureau et 
al. 1999; Rand & Wallin 2004; Hernandez et al. 2005; Emsellem 
et al. 2006; Fathi et al. 2007; Beckman et al. 2008; Meidt et al. 
2008), though the conditions under which it can be applied to 
the gas component are not yet fully understood. 

In this letter, we estimate the bar pattern speed of UGC 
628. This galaxy was selected because its morpological type 
and dark matter properties strongly differ from those of galax- 
ies usually observed for this measurement. Indeed, UGC 628 is 
a Magellanic barred, low surface brightness galaxy (LSB, here- 
after). Like most galaxies of this kind, its mass is entirely domi- 
nated by dark matter at all galactocentric radii (de Blok & Bosma 
2002). Our objective is to investigate how the dynamical prop- 
erties of a bar in a dark matter dominated galaxy compare with 
those of bars in more typical high surface brightness galaxies. 

For this measurement, we used optical 3D spectroscopic ob- 
servations of the ionized gas in UGC 628. It is a difficult task 
in the emission line gas and virtually impossible in the stellar 
absorption lines because the emission is very faint. This surely 
explains why pattern speeds have predominately been measured 
in high surface brightness galaxies. 

The kinematical and photometric observations are described 
in Section[T]and their analysis is done in Section[3] A discussion 
and brief concluding remarks are presented in SectionH] 
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2. Object Characteristics and Observations 

UGC 628 (D = 71.2 Mpc, / = 56°) was observed with 
FaNTOmM, an optical Fabry-Perot interferometer (Hernandez 
et al. 2003) mounted on the Canada-France-Hawaii Telescope 
(CFHT). These observations are part of a large survey of the 
kinematics of low surface brightness galaxies (see Chemin et al. 
in preparation for more details). The angular sampling of the Ha 
datacube is 0"49 per (square) pixel, corresponding to ~ 170 pc 
at the distance of the galaxy. The resolving power of the etalon is 
12450 at the redshifted Ha line. The spectral range of 262.7 km 
s"' was scanned through 24 channels. The exposure time was 
8.1 minutes per channel. The data reduction procedure is fully 
described in previous papers (e.g., Chemin et al. 2006). Figure[T] 
shows a 7?-band image of UGC 628 (de Jong 1996), the distri- 
bution and the velocity field of the Ha emission line. The broad 
band image shows a seemingly undisturbed spiral galaxy with 
an inner stellar bar aligned approximately North-South. The Ha 
emission displays a faint disc of diffuse gas. It is brighter in the 
bar region as well as along the two spiral arms which are more 
or less parallel to the major axis of the disc. 



3. Analysis 

3.1. The bar semi-major axis 

The intrinsic radius of the bar, aB, can be estimated following 
several ways (Athanassoula & Misiriotis 2002 and references 
therein). We use the two methods presented in Wozniak et al. 
(1995) and Aguerri et al. (2000) to derive ae for UGC 628. 

In the first method, isophotes are fitted to images of the 
galaxy. We used the BVRI images from de Jong (1996) and 
derived profiles of ellipticity and position angle, P. A., of the 
isophotal major axis using the ellipse task in iraf. The P.A. 
is defined as the angle between the semi-major axis of the re- 
ceding half of the velocity field and North (angle calculated 
North through East). Foreground and clumpy emission regions 
are masked before fitting. All photometric bands give similar re- 
sults though the B- and /-bands show more scatter which could 
be due to the low signal-to-noise ratio of the images. Results 
given hereafter have been obtained by averaging the values from 
all photometric bands. Only the results for the 7? -band image are 
shown for reasons of clarity. Wozniak et al. (1995) showed that 
ae can be estimated by determining the location of a discontinu- 
ity in the distribution of P.A. and ellipticity with radial distance. 
The change is such the regions in the bar show constant P.A. 



which then abruptly changes to the disc P.A.. Similarly, there is 
an abrupt change in the ellipticity from a fairly flat distribution to 
something significantly more round. Figure|2]shows the R-hmd 
radial profiles of P.A. and ellipticity of UGC 628 in the central 
regions of the galaxy. UGC 628 has a very small bulge (Fig.[T]) 
that should not contaminate the profiles significantly. The ellip- 
ticity peaks around a projected radius r = 9"4 before the end 
of the bar. The bar position angle is constant on the sky plane 
(~ 163°) and well separated from the disc orientation (P.A. = 
138°). The projected bar radius is irf3, which corresponds to 
13"2 after deprojection. 

A second method uses a Fourier decomposition of an image. 
Aguerri et al. (2000) define the bar as having the radial range 
where Ib/Iib > 0.5[max(lB/IiB) - min(lB/IiB)] + min(lB/IiB)- 
In this expression Ib = lo + I2 + I4 + l6 is the bar intensity. 
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Fig. 2. Results of isophotal ellipse fitting of the 7? -band image of 
UGC 628. The position angle and ellipticity of the isophotes are 
shown on top and bottom panels respectively as a function on the 
projected radius. The dashed vertical line indicates the projected 
bar semi-major axis. The shaded area represents the kinematic 
major axis of the disc at P.A. = (139 + 3)°. The horizontal line 
in bottom panel is the ellipticity of the disc for an inclination of 
56°. 
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with In, being the even Fourier coefficients of order m, and Iib = 
lo - I2 + I4 - le is the interbar intensity. We fitted the Fourier 
coefficients with a least-squares routine to a deprojected version 
of the 7? -band image of UGC 628. The results of the Fourier 
analysis is shown in Fig. [3] Following the prescription of Aguerri 
et al. (2000), the bar region extends between 2"5 < R < 14"8. 
The bar radius is the upper limit of this range. 

The two bar radius estimates are in very good agreement. 
The final bar semi-major axis is derived as the average of these 
two values, which gives ae - (14.0 + 0.8)" or (4.8 + 0.3) kpc. 

3.2. The bar pattern speed 

We use the Ha velocity field of UGC 628 (Fig. [Hand © with 
which intensity weighted mean positions and velocities are mea- 
sured along many slices aligned parallel to the kinematical major 
axis of the galaxy. This is the TW method. 

Though reasonable values of pattern speeds are derived us- 
ing Ha observations (e.g., Beckman et al. 2008), the conditions 
under which the method is applicable to the warm ionized gas 
is not entirely clear However, Hernandez et al. (2005), through 
the use N-body/SPH simulations found that the method was ro- 
bust for the ionized gas but with proviso that prominent regions 
where strong gas shocks occur should be avoided. Such regions 
introduce mean local velocity and position discontinuities that 
can bias slope measurements used in TW method. Fortunately, 
such discontinuites are easily identifiable in the TW diagrams. 

The kinematical position angle of the disc (P. A. = (139+3)°) 
defined the orientation of the slices. It is derived from a tilted- 
ring model of the Ha velocity field (Chemin et al., in prep.). 
There is very good agreement between the photometric and kine- 
matical major axes 03.1l l. The pattern speed of the bar is taken 
to be constant as a function of radius so that a single value can 
be derived using the TW equation of © We refer to Merrifield, 
Rand & Meidt (2006) and Meidt et al. (2008) for discussions of 
the radial dependence of pattern speed measurements in galax- 
ies. 

In practice, the integrals of Ev and 2x are not done within 
] - 00, +oo[ because of the finite extent of the galaxy due to the 
surface brightness detection Umit. The integrals are calculated 
along slices of maximal aperture X,nax, with the value of X,nax 
differing from one slice to another one. Results for two aper- 
tures values Xapeitui-e = X^ax - f and Xapeitui-e = X^ax are shown 
in Figure m where e is a small fraction of X^ax- Most of the in- 
tegrals have converged (or nearly so) to stable values <X> and 
<V> at these two apertures (filled symbols and solid lines). The 
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Fig. 3. Fourier analysis of a 7J-band image of the LSB galaxy 
UGC 628. The radial profile of the bar intensity Ib to the in- 
terbar intensity Iib is shown. The vertical dashed lines indicate 
the bar region. The outer line at R = 14.8" corresponds to the 
deprojected bar semi-major axis. 
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Fig. 4. Tremaine-Weinberg integrals as a function of the aperture 
length (or extent of the domain of integration). Filled (open, re- 
spectively) symbols and solid (dotted) lines are for the integrals 
that have (have not) converged to stable values. 
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Fig. 5. Bar pattern speed measurement of the LSB galaxy UGC 
628. Filled (open) symbols are for data inside (outside, respec- 
tively) the bar region. The solid line is the best fitting result to the 
filled symbols, with slope Op sin(i). The dashed line is the best 
fitting result to the open symbols, with a slope corresponding to 
the pattern speed of the spiral arms of UGC 628. 



TW method thus seems to be applicable to the ionized gas data of 
UGC 628. Integrals that have not converged within X,nax are dis- 
played with open symbols and dotted lines. They correspond to 
the edges of the gas distribution and velocity field, i.e. the most 
distant regions to the major axis. This is not surprising because 
the slices in these regions contain very few points. 
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Figure |5] shows the resuh of applying the Tremaine- 
Weinberg method to UGC 628. Only integrals which have con- 
verged within Xaperture = X^ax are drawn and kept in the slope 
fitting of the relationship. The uncertainties in the velocities are 
typically of the order of 3-4 km s ' . They depend mainly on the 
accuracy of the wavelength calibration of the raw datacube. The 
accuracy on the velocity centroid of an emission line in each 
pixel is indeed a fraction of a velocity channel, as derived from 
observations of a reference Neon lamp. 

Due to the lopsided galaxy morphology, the average posi- 
tions are preferentially observed at positive offests. A similar ef- 
fect has already been observed in other Ha or CO observations 
(e.g.. Rand & Wallin 2004; Hernandez et al. 2005; Emsellem 
et al. 2006). As a check for consistency with the stellar tracer, 
<X> values have been derived from the 7? -band image. They 
have been observed to be offset towards positive values as well. 

Only points which location falls within the bar region R < 
14" in the 2D Ha image are used to measure the bar pattern 
speed (filled symbols). The best fit gives Q.p = 11.3 km s"' 
kpc"^ This value does not change when a larger region is used, 
R < 14"8 and only increases by 2% when a smaller region is 
used, R < 13''2. All other points lying beyond the bar extent 
(open symbols) are considered to trace another pattern speed 
(e.g., that of the spiral arms in the disc). An amplitude of 8.8 
km s"' kpc~' is estimated for them, which value is lower than 
the bar pattern speed. 

The total uncertainty on Op is obtained by adding in quadra- 
ture the uncertainties in Qp due to the uncertainty in the position 
angle (3°), the inclination (2°), the extent of the region used for 
fitting Q.p (which is the error on the region of influence of the bar, 
i.e., fixed at the bar radius uncertainty of 0''8) and to the formal 
Icr uncertainty in the fit. All of the three later uncertainties are 
of-order 0.1 - 0.2 km s"' kpc"', which is negligible compared 
to uncertainty due to the uncertainty in the position angle. 3D 
spectroscopy allows us to derive Op at different position angles 
(between 136° and 142°) and we found a significant change of 
~18% in Qp (~2.0 km s"' kpc"'). Qp increases (decreases) to- 
wards smaller (larger) position angles than P. A. = 139°. This 
is similar to the effect seen by other authors (Debattista 2003; 
Debattista & Williams 2004). Because of this dominant source 
of uncertainty, we adopted Qp - (1 1.3 + 2.0) km s"' kpc"' as the 
likely best value. 



4. Discussion and conclusion 

Figure |6] (top panel) shows the Ha rotation curve of UGC 628 
(Chemin et al. in prep.), derived using a tilted-ring model of the 
velocity field. The angular velocity profile Q(R) = (bottom 
panel) was used to determine the corotation radius of the bar 
with the disc (Q = Qp). The corotation radius is Rc = 9. 8^2 
kpc implying a ratio - Rc/a^ - 2.0^^ 3. Given that observed 
bars having l.O <'R< 1 .4 are defined as being fast (e.g., Gerssen 
et al. 1999) and larger values of correspond to slow bars, UGC 
628 appears to host a slowly rotating bar. 

The slowest bar found before UGC 628 was for the blue 
compact dwarf galaxy NGC 2915 having H > 1.7 (Bureau et 
al. 1999). What is surprising here is that the bar-like structure of 
NGC 2915 is only seen in its neutral hydrogen distribution, not 
in the distribution of its stellar light. Purely gaseous bars are very 
rare in galaxies because bars are generally devoid of gas. Apart 
from that noteable exception, fast bars are always seen in galax- 
ies (Gerssen et al. 1999, 2003; Aguem et al. 2003; Debattista 



& Williams 2004; Corsini et al. 2007). To our knowledge, UGC 
628 is thus the first case of an obvious slow stellar bar in a galaxy. 

A striking point is that all fast bars are found in normal high 
surface brightness early-type discs. The stellar mass dominates 
the dark matter mass in their central regions, their disc is thought 
to be maximum. It is known that low surface brightness discs 
are dark matter dominated at all galactocentric distances, even 
with realistic mass-to-light ratio for their stellar disc (de Blok 
& Bosma 2002). From that, it seems that minimum disc galax- 
ies tend to host slow bars, in opposition to what is observed in 
maximum disc galaxies. Moreover, the fact that UGC 628 is at 
the opposite side of the Hubble diagram to the bright galaxies is 
perhaps the reason why it is observed a slow pattern speed for 
it. With two cases of slow stellar and gaseous bars in UGC 628 
and NGC 2915 (respectively), it is thus tempting to argue that 
the amplitude of bar pattern speeds decreases as a function of 
the Hubble type. This is in agreement with results from N-body 
simulations (Combes & Elmegreen 1993). 

Other numerical simulations explain how fast bars can be 
formed and maintained in a high surface brightness galaxy, pref- 
erentially within a maximum disc (Debattista & Sellwood 2000; 
Valenzuela & Klypin 2003). These simulations are in agree- 
ment with many observations of pattern speeds. Debattista & 
Sellwood (2000) argue that barred galaxies with moderately 
dense halos (like LSBs) should host a slow bar. Our result is 
therefore in full agreement with their expectation. Other numer- 
ical works investigate the bar evolution in LSB discs embedded 
in cosmological NFW halos (Mayer & Wadsley 2004). However 
they do not address the question of the bar pattern speed and the 
bar length compared with the corotation radius. 

Additional works are needed to confirm this new result. First, 
additional observations of barred low surface brightness discs 
would help to conclude whether UGC 628 is an isolated case 
or whether slow bars are common in this type of galaxies. They 
would also help to investigate whether pattern speeds really de- 
crease as a function of the morphological type. High-resolution 
Hi or optical interferometry like the current Fabry-Perot mea- 
surement seems to be the key for the Tremaine- Weinberg method 
because it is very difficult to observe stellar absorption lines 
of LSBs and the molecular gas is hardly detectable in them 
(O'Neil & Schinnerer 2004). A careful analysis of the datacubes 
of barred LSBs from our sample will be done. 

Finally, additional numerical simulations would also be very 
helpful to investigate the properties of bars in low surface den- 
sity, dark matter dominated galaxies. Low mass surface density 
discs are responsible of the cusp-core controversy (de Blok & 
Bosma 2002; Hayashi, Navarro & Springel 2007). Kinematical 
observations of LSBs indeed show that their dark halo seem to 
have a constant density in their core while steeper density pro- 
files are expected in Cold Dark Matter numerical simulations 
(e.g. Navarro, Frenk & White 1997). How would the pattern 
speed of a bar embedded within a cusp compare with that of 
a bar embedded within a core halo? Such simulations could per- 
haps help in constraining the shape of the inner density profile 
and alleviate the cusp-core controversy. 
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Fig. 6. Rotation curve (top panel) and angular velocity curve 
(bottom panel) of the LSB galaxy UGC 628. Dashed and dot- 
ted horizontal lines are the measured bar pattern speed of Qp - 
(11.3 + 2.0) km s"' kpc"'. Dashed and dotted vertical lines re- 
spectively indicate the bar semi-major axis a^, the corresponding 
corotation radius Rc and its lower and upper limits (Rein. Rcout) 
allowed by the uncertainties on Op. 
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